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The type Iln supernova 1994W: evidence for the explosive 
ejection of a circumstellar envelope 



3 



^1- 

o 



Nikolai N. Chugai 1 *, Sergei I. Blinnikov 2 , Robert J. Cumming 3 , Peter Lundqvist 
Angela Bragaglia 4 , Alexei V. Filippenko 5 , Douglas C. Leonard 6 , Thomas Matheson 5 ' 7 , 



^ ; Jesper Sollerman 3 



, 1 Institute of Astronomy, RAS, Pyatnitskaya 48, 109017 Moscow, Russia 
^■y 2 Institute for Theoretical and Experimental Physics, 117218 Moscow, Russia 

' ^Stockholm Observatory, Department of Astronomy, Stockholm University, AlbaNova University Center, SE-106 91 Stockholm, Sweden 
4 Osservatorio Astronomico di Bologna, via Ranzani 1, 40127 Bologna, Italy 
5 Department of Astronomy, University of California, Berkeley, CA 94720-3411 USA 
• 6 Five College Astronomy Department, University of Massachusetts, Amherst, MA 01003-9305 USA 
7 Harvard- Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 USA 

J> ' Accepted 12 May 2004. Received 11 May 2004; in original form 29 March 2004. 

ON ' 



ABSTRACT 



CO 

in 

We present and analyse spectra of the Type Iln supernova 1994W obtained be- 
\ tween 18 and 203 days after explosion. During the luminous phase (first 100 days) 

— -» i the line profiles are composed of three major components: (i) narrow P-Cygni lines 

*-£h ' with the absorption minima at —700 km s -1 ; (ii) broad emission lines with blue ve- 

locity at zero intensity ~ 4000 km s -1 ; and (iii) broad, smooth wings extending out 
Q ■ to at least ~ 5000 km s -1 , most apparent in Ha. These components are identified 

S-^ | with an expanding circumstellar (CS) envelope, shocked cool gas in the forward post- 

r/2 . shock region, and multiple Thomson scattering in the CS envelope, respectively. The 

absence of broad P-Cygni lines from the supernova (SN) is the result of the forma- 
tion of an optically thick, cool, dense shell at the interface of the ejecta and the CS 
envelope. Models of the SN deceleration and Thomson scattering wings are used to 
y\ ' recover the density (n « 10 9 cm -3 ), radial extent [~ (4 — 5) x 10 15 cm] and Thomson 

optical depth (tt 2.5) of the CS envelope during the first month. The plateau-like 
SN light curve is reproduced by a hydrodynamical model and is found to be powered 
by a combination of internal energy leakage after the explosion of an extended pre- 
supernova (~ 10 15 cm) and subsequent luminosity from circumstellar interaction. The 
pre-explosion kinematics of the CS envelope is recovered, and is close to homologous 
expansion with outer velocity ~ 1100 km s -1 and a kinematic age of ~ 1.5 yr. The 
high mass (~ 0.4 Mq) and kinetic energy (~ 2 x 10 48 erg) of the CS envelope, com- 
bined with small age, strongly suggest that the CS envelope was explosively ejected 
~ 1.5 yr prior to the SN explosion. 

Key words: supernovae - circumstellar matter - stars: supernovae: individual (SN 
1994W) 



1 INTRODUCTION 

Recent studies of supernovae (SNe) have made considerable 
progress in our knowledge of what makes a star explode. 
From observations of light curves and spectra, theories of 
pre-supernova evolution are now better constrained than 



ever. Conventional wisdom says that Type II SNe (SNe II) 
are caused by core collapse in massive, usually red super- 
giant stars. In general, theory has no trouble accounting for 
the main features of the spectra and light curves of these ob- 
jects. However, some SNe II are remarkably different, and 
many of their observational properties are not yet under- 
stood. 
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These events, sometimes known as narrow-lin e Type 
II SNe or SNe Iln JSchlegellll99ol: lFilippenkolll997lh show 
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in early spectra the presence of strong, narrow Balmer 
emission lines on top of broad emission l ines. Exam- 
ples of prominent SNe Iln include 1978K llRvder et alJ 
Il993l) SNe 1983K iNiemela. Ruiz fc Phillipd [19851) . 1984E 



JPopita et al]ll984L 1986J llRupen et al.lll987tlLeibundgut 



1991) 



1994) 



1987F 
1988Z 



Filippenk 
Fflrpjsenki 



1991; Turatto et al 



1993; 




IWegner fc Swanson 

Stathakis fc S adler 
' I I ■ 1 



iBenetti et al 



1998), 



1994W 



■Cumming. Lundavist fc Meikld 1 1994 [Meikle et alJ Il994t 
ISoll^nriM i^Qrmrnhi^^ ^un£krvislJll998l , he reafter SCL98), 



1995G JPastorello et alJl2002T> and 1995N jFransson et al] 



are SNe 1997cv iGermanv et al.lll994l: iTuratto et al.ll200(J). 


1998S (Filippcnko & Moran 1998; Bowen et al. 


2000; 


Leonard et alJ 


2000: iFassia et alj|200ll: IChugai et alJ 


2002; 


Fransson et al. 


2004), 1999E (Filippcnko 2000; Rigon ct al. 


2003) and, most remarkably, SN 2002ic, which has been 



identified as a likely SN la in a dense circumstellar (CS) 
envelope iHamuv et al.ll2003j) . are more recent examples of 
SNe with narrow lines. 

The general wisdom is that the narrow lines of 
SNe Iln originate from the ionized, dense circumstellar gas 
jHenrv fc Branch! Il987t lFilippenkolll99lll . The interaction 
of SN ejecta with the dense CS gas modifies the SN op- 
tical sp_ectouinandjnay fully power the SN Iln luminos- 
ity <IChugailll990l 1992). Interaction with the dense CS gas 
is indicated also by strong radio and X -ray flux detecte d 
in some SNe Iln. SNe Iln are diverse jFirippenkolliggfl) . 
probably reflecting variations in CS gas density and struc- 
ture (smooth vs. clumpy) and SN ejecta parameters (mass 
and energy). On top of this there is also the possibility of 
asymmetry of both the ejecta and CS gas. Nevertheless, 
the primary factors responsible for the diversity (mass of 
main-sequence star, explosion mechanism, structure of pre- 
supernova, mass-loss mechanism and history, and progenitor 
binarity) are still unknown. All these uncertainties, along 
with the possibility of using them to probe pre-supernova 
behaviour prior to core collapse, add to the interest in these 
phenomena. 

SN 1994W, discovered on 1994 July 29 (UT dates 
are us ed throughout this paper) at the pre-maximum 
phase JCortini fc Villilll994l) . is among the brightest known 
SNe Iln. Its proximity (25±4 Mpc; SCL98) and early dis- 
covery made it an ideal case for detailed study. Preliminary 
interpretation of the spectra showed that SN 1994W ex- 
ploded in the dense CS envelope with a characteristic ra- 
dius of ~ 10 15 cm, while the amount of ejected 56 Ni was low 
(< 0.015 M Q ; SCL98). 

Here we further advance our understanding of both the 
CS envelope properties and the phenomenon of SN 1994W as 
a whole by analysing the spectra and photometry discussed 
by SCL98 as well as other, hitherto unpublished, spectra. In 
our study we take two different approaches to modelling the 
data. The first is the simulation of the Ha profile, which pro- 
vides an efficient and straightforward probe of the Thomson 
optical depth and density of the CS envelope. The second ap- 
proach is a hydrodynamical simulation of the SN explosion 
and the light curve based upon the upgra ded version of the 
code STELLA ( Blinnikov ct al 1 ll998Ll200(tl . An early version 
of this code proved effective in computing the hydrodynam- 
ics and light curve for int eraction with a dense CS env elope 
in the case of SN 1979C jBlinnikov fc Bartunovlll993h . 



This paper has the following structure. We start with 
a description of the data (Sections [5] and [3J and the gen- 
eral picture which arises from a qualitative analysis (Section 
Next we quantify parameters of the CS envelope using 
a thin-shell deceleration model and a line-profile simulation 
(Section|5Jl. We then use the gross parameters of the recov- 
ered density distribution as input for detailed hydrodynam- 
ical modelling ( Section [SJ. Our findings are summarized and 
discussed in Section Q 



2 OBSERVATIONS 
2.1 Photometry 

The photometric data were compiled and described in 
SCL98. In our light-curve models in Sectional we use the 
photometry up to day 197, assuming SN 1994W exploded 
on 1994 July 14.0 (SCL98). For a discussion of the explosion 
date, see Section [3] 



2.2 Spectroscopy 

The spectroscopic observations are described in Table 
Some of these, taken from La Palma, were presented in Ta- 
ble 2 of SCL98. We supplement them with spectra taken at 
the Lick, Keck, and Bologna Observatories. The data cover 
epochs 18-203 days after explosion. We now describe these 
spectra in detail. 

The La Palm a spectra were reduced using the figaro 
JShortridgd Il99fj) and NOAO iraf 1 packages. The CCD 
frames were bias-subtracted, flat-fielded, and corrected for 
distortion and cosmic rays. Wavelength calibration was car- 
ried out using arc spectra of copper-neon and copper-argon 
lamps, and should be accurate to better than 0.1 A for the 
high-resolution WHT and INT spectra, 0.5 A for the low- 
resolution WHT and INT spectra, and 1 A for the NOT 
spectrum. 

The Bologna spectrum was obtained with the Bologna 
Astronomical Observatory 1.5-m telescope with the spec- 
trograph BFOSC on 1994 July 31.9 and reported by 
iBragaglia et~aH dl994l) . 

Spectra were also obtained at the Shane 3 m reflec- 
tor at Lick Observatory on days 21, 49, 57, 79 and 89, and 
at the W. M. Keck II 10 m telescope on day 197. The 
Lick obse rvations were made w ith the Kast double spec- 
trograph jMiller fc Stond 119931) . The Keck spectrum was 
obtaine d using the Low Resolution Imaging Spectrometer 
(LRIS; (lOke et al.lll993l) ) with a 1200 lines mm" 1 grating. 
For details, see Tabled 

For the Lick data, all one-dimensional sky-subtracted 
spectra were extracted in the usual manner, and each spec- 
trum was then wavelength and flux calibrated, as well as 
corrected for continuum atmosp heric extinction and tel- 
luric absorption band s llWade fc Horndll988l : lBesselllll999l : 
iMatheson et al.ll2000h . 



1 IRAF (Image Reduction and Analysis Facility) is distributed by 
the National Optical Astronomy Observatories, which are oper- 
ated by the Association of Universities for Research in Astronomy, 
Inc., under cooperative agreement with the US National Science 
Foundation. 
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Table 1. Log of spectroscopic observations 



Date 








Phased 


Telescope / 


Spectral 


Flux 


Seeing 


Wavelength 


Observers 


CUT) 










spectrograph 


res. (A) 


standard 


(arcsec) 


coverage (A) 




1994 


07 


31, 


9 


17.9 


RAO /RFOSO 


13 


Rr>4-25°3941 


1.6 


4600-7800 


AB 


1994 


07 


31. 


9 


17.9 


INT/IDS 


0.25 




2.6-3.0 


6507-6675 


EZ 


1994 


08 


01 




21.5 


Tjirk /Kast 


7 


RD-LI 7° 47ns 


2.5 


4250-7020 


AIR AF CP 


1994 


08 


13. 


9 


30.9 


WHT /TSTS 


2.7 


SP1 337+705 

ijl -LOO 1 \^ I 


1.1-1.7 


4290-8005 


MB 


1994 


09 


01 




49.5 


T.irV /K"ac;t 

J_J 1 K,iS. 1 Iv do L 


5 


RD+26°2606 


44- 


qi 90—1 runn 


A TR AF TM 

rtij _L> ; . \ I . _L 1V±. 


1994 


09 


08. 


9 


56.9 


WHT/ISIS 


11 


Grw+70°5824 


0.9-2.9 


3300-9360 


RR 


1994 


09 


08. 


9 


56.9 


WHT/ISIS 


2.7 


Grw+70°5824 


1.0-1.5 


4255-5060, 
6345-7150 


RR 


1994 


09 


09 




57.5 


Lick/Kast 


4 


RD+28°4211, 
RD+26°2606, 
RD+17°4711 


3-5 


3140-9920 


AJB, AF, TM 


1994 


10 


01 




79.5 


Lick/Kast 


4 


Feige 110, 
RD+26°2606 


2-3 


3130-8040 


AF, LH, TM 


1994 


10 


11 




89.5 


Lick/Kast 


4 


Feige 110, 
BD+26°2606 


3-4 


3160-8020 


AF, TM 


1994 


11 


12. 


3 


121.3 


INT/IDS 


6 




1.6 


5485-7370 


RR, NO 


1995 


01 


26. 


6 


196.6 


Keck/LRIS 


6 


HD 19445 


1.2 


5628-6939 


AF, LH 


1995 


02 


01. 


9 


202.9 


NOT/LDS 


12 


Feige 34 


0.9 


5000-10000 


JS, RC 



Note: (a) Days after 1994 July 14.0. (b) Observers: AB - A. Bragaglia; AF - A. Filippenko; AJB - A. J. Barth; CP - C. Peng; EZ - E. 
Zuiderwijk; JS - J. Sollerman; LH - L. Ho; MB - M. Breare; NO - N. O'Mahony; TM - T. Matheson; RC - R. Cumming; RR - R. 

Rutten. 



The spectra were flux-calibrated by comparison with 
the standard stars listed in Table Q The sky position of 
SN 1994W in late 1994 meant that the observations were 
made when the supernova was at rather low altitude, and 
correspondingly high airmass. We corrected for the differ- 
ence in atmospheric extinction brought about by the differ- 
ing zenith distances of the supernova and the standard star. 
Airmasses at diff erent zenith dis tances were taken from the 
relation given by iMurravl il983Jh and we used the tables of 
ex tinction coef ficient versus wavelength for La Palma given 
bv lKinsNl985th 

With two exc eptions, all our sp ectra were taken at the 
parallactic angle l|FiliDPenkol 11982) . minimising slit losses 
due to atmospheric dispersion. The first exception was the 
day 18 Bologna spectrum of SN 1994W, which shows con- 
siderable loss of flux in the blue. We corrected the slope of 
the spectrum to match the B and V photome try from the 
same night re ported bv i Bra gagliaet aljlll994h using filter 
functions from lOualandi fc Merighil l|200lf) . 

The second exception was the standard-star observation 
for the day 31 La Palma spectrum. While the SN 1994W 
spectrum was taken at the parallactic angle, the standard- 
star spectrum was taken at an airmass of 1.6, with the slit at 
position angle 33° when the parallactic angle was 119°. No 
contemporary filter photometry is available for correcting 
the supernova spectrum. The effect on the day 31 supernova 
spectrum depends on how the standard star was positioned 
in the slit, which in turn depends on the combination of the 
spectral response of the acquisition camera at the telescope. 
The acquisition camera at the WHT in 1994 August was 
a Westinghouse ETV-1625 whose response curve peaked at 
4400 A, with half-power points at 3400 A and 6200 A (C. 
Jackman & D. Lennon, 2004, private communication). The 
standard star, Grw +70° 5824, is a DA3 white dwarf with 



a blue spectrum. It therefore seems likely that the slit was 
positioned so that maximum transmission was at around 
4400 A, at the blue end of the spectrum. 

We have estimated the amount of flux that would be 
lost with wavelength, and find that the losses would largely 
have been at the red end of the standard spectrum — as 
much as 80 per cent of the flux at 8000 A. The effect on 
the calibrated supernova spectrum was thus to make it ap- 
pear redder than it really was. Since the day 31 spectrum is 
already remarkably blue (Section 13.1. 311 . we have made no 
further correction to the spectrum; the calibration can be 
regarded as a conservative estimate of the already extreme 
continuum slope. 

No observations of flux standards were available for the 
INT spectra taken on days 18 and 121. A very rough flux- 
calibration of the day 18 INT spectrum was made by assum- 
ing that the flux in the wings of the Ha line (at the red and 
blue edges of the spectrum) declined at the same rate as the 
visual magnitude estimates between days 18 and 31, i.e., by 
a factor of about 1.6 in fx- A similarly rough calibration 
of the day 121 spectrum was made by comparing with our 
i?-band photometry from day 123 (SCL98). 

The spectrum taken on day 197 was scaled so that the 
flux in the narrow Ha feature matched that in the day 203 
spectrum. The continua then match well in the overlap re- 
gion. 

Spectra from days 18, 57, 121 and 203 were displayed in 
SCL98 (their Figs. 2 and 3). They also showed the portion 
of the spectrum on day 31 covering Ha. 

As a final correction, we noted that all the spectra cover 
the entire V band, and we scaled each spectrum to match 
contemporary V-band photometry. The procedure was as 
follows. We assumed that the day 57 La Palma spectrum was 
correct. Next, we interpolated the V magnitudes in the light 
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Figure 1. Spectra of SN 1994W: days 18-203. The spectra have been shifted vertically for clarity. No heliocentric correction has been 
applied to the wavelength scale. The ticks on the right-hand side mark the zero level for each spectrum. The spectra from day 79 and 
later have been multiplied by a constant, noted in parentheses. The La Palma spectra from days days 18 and 57 and the day 197 Keck 
spectrum are not shown. Telluric absorption lines in the day 31 spectrum are marked. 



curve presented by SCL98. This gave us relative fluxes in V 
for all the epochs. We then integrated all the plateau spec- 
tra (days 21-89), mu ltiplying them by a standard V filter 
function jAllenlll973l) with transmission values interpolated 
between points at 100-A intervals. The whole spectrum was 
then scaled linearly in flux so that the integrated V magni- 
tude corresponded to the V-band photometry in SCL98. 



the spectral development for each epoch in turn, and then 
concentrate on the individual features and ionic species. 

In this section, quoted line fluxes (see Table were 
measured by interpolating the neighbouring continuum 
across the line in question. The line fluxes are net fluxes 
and do not correct for the narrow P-Cygni absorption com- 
ponents. 



3 RESULTS 

During the period of our observations, the spectrum of SN 
1994W (Figs. Q and |2J was dominated by narrow Balmer 
emission lines, initially with P-Cygni absorption features, 
and accompanied by a wide variety of weaker lines show- 
ing broad emission and /or narrow P-Cygni profiles (cf. 
iFilippenko fc Bartbll9"94l) . Line identifications are presented 
in Tableland in Fig. [3] In the following sections, we examine 



3.1 Spectral evolution 

3.1.1 Day 18 

The low-resolution Bologna spectrum showed a blue contin- 
uum with prominent lines of Ha and H/3, which are nar- 
row with broad wings. The latter have been interpreted as 
an effect of multi ple Thomson s cattering of line photons in 
the CS envelope (IChugaill200"lh . Our high-resolution spec- 
trum for this epoch (Fig. 3 in SCL98) covers only the Ha 
line, but shows a narrow P-Cygni absorption at about —800 
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Figure 2. Spectra of SN 1994W: days 18—121. Same as Fig.Q but with the ordinate stretched to show the weaker features more clearly. 
The day 197 and 203 spectra are omitted. Telluric absorption lines in the day 31 spectrum are marked. 



km s _1 not visible in the low-resolution spectrum. The emis- 
sion peak is close to the velocity of nearby H II emission, 
measured at 1249±3 km s _1 , which we adopt as our best es- 
timate of the heliocentric velocity of the supernova. The in- 
terstellar Na I D lines give a somewhat lower value, 1185T15 
km s~\ measured from the day 31 spectrum. 



3.1.2 Day 21 

The rapid brightening of the supernova at this epoch 
(SCL98) was reflected in increasingly strong Balmer lines. 
The spectrum showed a strong blue continuum, and promi- 
nent lines of Ha, H/3 and H7, all exhibiting narrow P- 
Cygni profiles with broad emission wings. The Balmer decre- 
ment was remarkably shallow, with flux ratio Ha:H/3:H7 
= 1.1:1:0.5 (dereddened assum ing Eb-v =0.17 mag from 
SCL9 8 and the reddening law of lCardelli. Clayton fe Mathisl 
1989). The Balmer lines were accompanied by an asymmet- 
ric triangular broad emission feature corresponding to He I 



A5876, possibly blended with the Na I D doublet, and many 
narrow P-Cygni lines of Fe n. 



3.1.3 Day 31 

By day 31, the blue continuum appeared to have steepened. 
The calibration of the spectrum is uncertain, as noted in 
Section |2] above, but the evidence suggests that the spec- 
trum really was as remarkably blue as it appears to be. The 
Balmer lines still show a flat decrement, as on day 21, with 
measured ratios 1.2:1:0.6. The blue wings in H/3 and H7 ap- 
pear to be stronger relative to the red wings, compared to 
day 21 (Fig. [3] right-hand panel). 

Broad, triangular He I emission lines with full- width at 
half-maximum velocity width («fwhm) around 2500 km s _1 
were stronger than on day 21 (Figs. [2] and 0. He I A5876 
was accompanied by strong He I A7065. Broad He I A6678 
appeared to fill out the red wing of Ha, and the triplet 
line at A4471 is present as broad emission which underlies 
neighbouring Fe n lines. All but He I A7065 also showed 
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Figure 3. Line identifications for days 49 (upper spectrum) and 79. The spectra have been continuum-subtracted and shifted for clarity. 
Each line is indicated by a broken vertical bar marking heliocentric velocities —700 km s — 1 and km s — 1 . Horizontal bars longward of 
6700 A mark regions where telluric absorption has been removed. 



weak, narrow P-Cygni features with minima around —700 
km s _1 . 

A similarly triangular emission feature is seen at 7890 A, 
which we identify with the Mg II triplet at 7877-7896 A (see 
Section I3.3.4H . 

Stronger narrow P-Cygni lines with absorptions around 
-600 km s _1 are seen in O I A7773, Si II AA6347, 6371, and 
a large number of Fe II lines. The O I and Si II lines are dom- 
inated by absorption, while in Fe II the equivalent widths of 
the absorption and emission components are comparable. 



3.1.5 Day 57 

Eight days later, the Balmer decrement showed no signifi- 
cant change at 1.5:1:0.6. The spectrum showed ever stronger 
Fe II lines. No unambiguous trace of He I remained (Fig.yJ, 
suggesting that the emission feature at ~6000 A was now 
only due to Na I D. The emission in this feature broadened 
and flattened out, extending from —4000 km s _1 to +6000 
km s _1 (Figs. H and EJ. The Mg II A 7877-7896 feature was 
still present, but appeared more flat-topped than on day 
49. P-Cygni lines of O I A8446 and the Ca II near-infrared 
triplet, also weakly present on day 49, were clearly detected 
in the red. 



3.14 Day 49 

The day 49 spectrum (Figs.Q|21and|3J shows a redder con- 
tinuum. The Balmer decrement steepened to 1.5:1:0.6. The 
He I A5876 and Na I D blend was again asymmetric as on 
day 21, and appeared to have a broader blue wing. The only 
other remaining He I line, A7065, is no longer as prominent 
as on day 31. The Mg II A7877-7896 feature was still clearly 
visible. Similar broad features were detected at ~8230 A and 
~9220 A, lending support to the identification with Mg II. 

Many new Fe II lines are apparent, plus Ca II H & K 
and possibly the Ti II multiplet at around 3500 A. The wider 
wavelength coverage picked up a large number of overlapping 
narrow P-Cygni lines in the blue; most of the identified ones 
are H I and Fe II. 



3.1.6 Day 79 

The continuum appeared again to have reddened (Fig. 0. 
The Balmer lines decreased in strength relative to the now 
very numerous Fe II P-Cygni lines (Figs. |2] and The 
Ha/H/3 ratio was now 2.7±0.3, close to Case B. Narrow over- 
lapping absorption and P-Cygni features completely domi- 
nated in the blue. In the red, the feature at 7890 A had 
disappeared, but O I A8446 and the Ca II triplet had in- 
creased in strength. The narrow absorption in Na I D had 
increased markedly in strength since day 57, and broad Na I 
emission was no longer apparent. 

Narrow P-Cygni profiles of O I A7002 and Sc II ap- 
peared, the latter possibly accompanied by a broad emission 
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Table 2. Line identifications for days 18-121. Symbols are as follows: a — line seen in absorption only; b — broad emission, e — 
emission; p — P-Cygni pr ofile; — (dash) — no feature de tecte d, : (colon) — uncertain identification. The Fe II lines were identified with 
the help of the line li sts of lFuhr. Martin fc Wiesel ll988l) and lSigut fc Pradhanl 120031) . Wavelengths in A are from these references and 
from Ivan Hooa 11999). 
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is obvious that a dramatic change had occurred since day 
89, paralleling the supernova's precipitous fading after day 
~110 (SCL98). 

As reported by SCL98, three narrow emission lines are 
unambiguously detected: Ha with a upwhm of 780 km s , 
plus much fainter Na I D (dfwhm = 775 km s _1 ) with a P- 
Cygni absorption component, and Ca 11] A7291 (wpwhm = 
475 km s" 1 ) close to the red edge of the spectrum. The emis- 
sion peak of Ha had moved farther to the red than on day 
89, at about +70±100 km s" 1 . The Sc 11 P-Cygni feature at 
5530 A appears also to have survived since day 89, showing 
a sharp absorption feature at —800 ± 100 km s _1 . 

Underlying the narrow features is what appears to be 
a low, undulating continuum. SCL98 argued that the sharp 
drop to zero flux at 5620 A was evidence that the contin- 



component. Semi-forbidden lines of Fe 11] AA5991, 6150 were 
clearly detected for the first time. 

3.1.7 Day 89 

The spectrum on day 89 was essentially the same as on day 
79, but was weaker and showed a redder continuum. The 
Ha/H/3 ratio increased to 5.2±1.3, and H7 was now blended 
with Fe 11 lines. Both Na I D and O I A7002 continued to 
increase in strength. 

3.1.8 Day 121 

Our spectrum for day 121 (Figs.Qand|2j has smaller wave- 
length coverage and rather poor signal-to-noise ratio, but it 
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Figure 4. Luminosity evolution of various lines. Line fluxes are 
detailed in Tabled A distance of 25.4 Mpc and an extinction of 
Eb—v = 0.17 mag were adopted. 
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Figure 5. Selected line profiles; day numbers are indicated. The 
left-hand panel shows examples of narrow P-Cygni profiles, the 
middle panel shows broad lines (note the extended wings in Na I 
on day 57), and the right-hand panel shows the wings of the H I 
lines on day 31. The dotted vertical lines mark velocities of and 
-700 km s -1 . 



uum was rather composed of broad, overlapping emission 
lines. The steepness of the drop at 5620 A indicates that the 
velocity widths of such lines can be no greater than ~5000 
km s _1 . Similar feature s have been seen in la t e-time spectra 
of th e SN Iln 1997cy iGermanv et alJll994t iTuratto et all 
2000). A sharp drop at this wavelength is in fact also typ- 
ical of SNe la and some SNe Ib/Ic at late times, when 
the emission is dominated by lines of Fe II and [Fe n] (see 



lAxelrodlll980l) . These similarities suggest that SN 1994W 
had reached its nebular phase at this epoch. 

3.1.9 Day 197 

Our spectrum, not shown in Figs. and |5J detected only 
Ha, which was perhaps just resolved. Its velocity width was 
vfwhm ~ 625 km s _1 , consistent with no change since day 
121. Na I emission at the same strength relative to Ha as 
on day 121 would not have been detected, given the low 
signal-to-noise ratio of the spectrum. 

3.1.10 Day 203 

The day 203 spectrum (Fig. confirms that narrow, un- 
resolved Ha is the only detectable feature from the su- 
pernova. The galaxy background probably dominates the 
remaining, noisy continuum, though photometry indicates 
that the supernova can account for up to 50 per cent of the 
flux (SCL98). SCL98 put limits on the flux in broad lines of 
width ~4500 km s" 1 . 

3.2 General spectral appearance 

3.2.1 Line profiles 

The spectra reveal two major types of line profile-forming 
components (SCL98; Fig. [^J: narrow P Cygni profiles and 
broad lines. However, careful inspection reveals that while 
the blue velocity at zero intensity (BVZI) of the broad lines 
does not exceed —5000 km s _1 , Ha at the luminous stage 
shows a red wing extending at least up to +7000 km s _1 . 
This wing is clearly detectable in Ha on days 31, 49 and 
57 (Figs. [5] and We believe that this wing is caused by a 
Thomson scattering effect in the expanding CS envelope (see 
also Section 2J. Since the intensity of the Thomson scatter- 
ing wing is proportional to the line intensity, the wing may 
be present but not obvious in other, weaker lines. 

The intrinsically broad lines typically have a triangular 
profile. On day 31 they show an apparent skewing toward 
the blue, possibly due to an occultation effect. The Na I 
doublet on day 57 shows a flat-topped profile with BVZI~ 
4000 km s _1 . No other line shows such a flat-topped profile, 
nor does Na I doublet at other epochs, though the feature 
is dominated by He I A5876 earlier than this. This suggests 
that the material emitting Na I has a narrower distribution 
of velocity than the other broad lines (see Section [4.211 . 

3.2.2 Line strengths 

Assuming a distance of 25.4 Mpc to NGC 4041 and E B -v = 
0.17 mag (SCL98), we have plotted the luminosity evolution 
of the strongest emission lines in Fig. [I] Note the quasi- 
exponential decline of the Ha luminosity with exponential 
lifetime ~ 21 days. This fast decline is in dramatic con- 
trast to other SNe Iln, for example SN 1995G, whose nar- 
row Ha luminosity drop ped by only a fact or of three during 
the first 600 days iPastorello et alJl2002l) . The fast decay 
of the Ha luminosity is consistent with a small amount of 
56 Ni in SN 1994W (SCL98) and also indicates the absence 
of dense CS gas at large distances from the pre-supernova 
(see Section 0. 
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Table 3. Flux measurements of emission lines, in units of 10 14 erg s 1 cm 



Day 



Line 


18 


21 


31 


49 


57 


79 


89 


121 


203 


H 7 4340.46 




30(2) 


47(4) 


34(1.5) 


24.5(2) 










H/3 4861.32 


30(10) 


61(3) 


91(5) 


60(5) 


47(2) 


10(1) 


4(1) 






Fe II 5018.44 




1.7(0.1) 


0.7(0.1) 


2.1(0.5) 


2.7(0.7) 


2.0(0.7) 


1.6(0.2) 






He I 5875.62 




















+Na i< a ) 5889-95 




13.6(2) 


14(4) 


12(1) 


5(1) 


<1 


<1 






Na i( b > 5889-95 










0.14(0.05) 


0.2(0.02) 


0.6(0.1) 


0.11(0.03) 


<0.03 


Ha 6562.80 


25(5) 


82(3) 


130(6) 


110(2) 


84(2) 


32(3) 


21(1) 


2.2(0.2) 


0.1(0.025) 


He I 7065.22 


<2 




9.7(3) 


2.9(0.7) 


<0.7 


<0.5 


<0.2 






Mg rrM 7877-7896 






10(2.5) 


5.2(0.5) 


5.0(0.3) 


<0.6 


<0.6 






Mg II 9218-44 










6.3(1) 


7.5(0.8) 


1.0(0.3) 







Notes: (a) blend, (b) emission component only. 
Error estimates are la. Upper limits are 3a. 



3. 2. 3 Continuum 

We carried out black-body fits to the continuum. We obtain 
colour temperatures of 13000±2000 K on day 21, 15000 K on 
day 31, 10000 K on days 49 and 57, 7200 K on day 79 and 
7200 K on day 89. We find that the black-body fits work 
best for Eb-v = 0.15 mag, consistent with our estimate 
from Na I absorption (0.17 ± 0.06 mag; SCL98). 



3.3 Elements identified 

3. 3. 1 Hydrogen 

We detect the Balmer series from Ha to at least Hll, and 
possibly also members of the Paschen series from Pa<5 to 
Pall (Fig. The profiles seem to be a combination of 
broad emission line with BVZI > 4000 km s" 1 and a narrow 
P-Cygni absorption component (~ 1000 km s _1 ). The ra- 
tio of narrow-to-broad emission component decreases from 
Ha toward the higher Balmer lines. This indicates a nearly 
normal Balmer decrement for the narrow component and 
a flat or possibly inverse Balmer decrement for the broad 
component (Fig. |KJ . 

A flat or inverted decrement suggests strong deviation 
from the recombination case, related to collisions and ra- 
diative transitions in the field of the photospheric radia- 
tion. The observed line ratios suggest an excitation tem- 
perature for hydrogen's excited levels of around 15000 K. 
Two different thermalization mechanisms for hydrogen lev- 
els may operate in this situation: the level populations may 
be controlled by radiative transitions in the field of the 
external black-body photospheric radiation, or thermaliza- 
tion may be primarily collisional. Both mechanisms have 
been exploited to account for flat or inv erted Balmer decre- 
ments observed in cataclysm ic variables IIElitzur et alJll983t 
IWilliams fe Shipmanl Il988l) . Both thermalization mecha- 
nisms are plausible in SN 1994W as well. 



3.3.2 Helium 

We observe broad He I lines with triangular emission profiles 
and narrow P-Cygni features. The He I lines disappear after 
day 49 (Figs. Handel Tabled . He I A5876 and A7065 appear 



only in emission. He I A4471 and A6678 show weak absorp- 
tion features. A high ratio /(7065)//(5876) « 0.5 indicates a 
large contribution of collisional and/or radiative excitation 
compared to the purely recombination case. Unfortunately, 
this fact cannot be used as a straightforward indicator of 
electron concentration in the line-em itting region. Computa - 
tions for an extended parameter set iAlmog fe Netzerlll989l) 
show that this particular ratio may be reached in a wide 



range of concentrations (n e 



10 4 - 10 14 cnT 



with pro- 



nounced dependence on the optical depth in He I A3889. 



3.3.3 Oxygen 

O I A7773 is present with an absorption-dominated P-Cygni 
profile from as early as day 18 to day 89. From day 79, 
we tentatively identify O I A7002 with a P-Cygni profile. 
This line has not to our knowledge been seen in a supernova 
spectrum before, b ut was identified in the proto-planetar y 
nebula Henize 401 llGarcia-Lario. Riera fe Manchaddll999T) . 
whose low-ionization spectrum bears some similarity to that 
of SN 1994W at these phases. 



3.3.4 Other metals 

Ca 11 H & K were strong in absorption at all epochs when 
we covered them. The near-infrared triplet lines were weak 
on day 49, but strengthened dramatically by day 79 (Fig. 01. 
This may reflect an increase in the Ca 11/Ca III ratio with 
the drop in radiation temperature at late epochs. 

The Na I D doublet's development is remarkable (Figs. 
0121 and |^J, even taking into account blending with inter- 
stellar absorption and in the earlier spectra with He I A5876. 
Up to day 49, a flat-topped profile emerged as He I faded. 
From day 57 on, the strength of the circumstellar P-Cygni 
absorption increased, presumably as a result of increasing 
Na 1 ionization fraction. 

Magnesium lines appear both in absorption and emis- 
sion. Narrow absorption with weak P-Cygni emission is 
seen in Mg 11 A4481 from day 21 to day 57. We iden- 
tify Mg 11 A7877-7896 as the source of the broad feature 
around 7890 A. As an alternative interpretation, we consid- 
ered [Fe xi] A7892. Strong [Fe xi] might be expected to be 
accompanied by emission in [Fe x] A6374, and there may 
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indeed be a broad feature underlying the Si II absorption 
lines to the blue of Ha. However, this feature does not per- 
sist to later epochs as the 7890 A feature does, nor are any 
other high-ionization lines identified at any epoch. In addi- 
tion, emission lines with similar widths to the 7890 A feature 
are seen at wavelengths corresponding to Mg II transitions 
around 9230 A, and possibly at 8200 A and 10900 A. Never- 
theless, we find it somewhat surprising that the broad emis- 
sion in Mg II AA7877-7896 is not accompanied by similar 
emission in O I A7773. 

The doublet of Si II AA6347, 6371 is present from days 21 
to 89, showing absorption-dominated P-Cygni profiles. The 
lines weakened relative to other features on days 79 and 89. 
Si II A6240 therefore seems to be an unlikely contributor to 
the feature at 6260 A (see below), which increases on these 
dates. 

Sc II em ission lines at around 5 530 A and 6250 A were 
identified by IPastorello et all i2002T) in the spectrum of the 
Type Iln SN 1995G. We identify both Sc II A5527 and A6246 
in our spectra, and Sc n may contribute to the feature at 
6260 A. In particular, A5527 shows some evidence of a broad 
emission component (Figs.|21and|3J. Its narrow P-Cygni core 
apparently persisted until day 121. 

We identify a host of Fe II lines (Fig. 01 Table all of 
which show narrow P-Cygni profiles. They may also be ac- 
companied by broad emission as in the case of other species, 
since the spectra show hints of underlying broad emission 
in regions of the spectrum where there are many Fe II lines 
(Fig. 0. By days 79 and 89, a few lines of forbidden and 
semi-forbidden Fe II are identified. These have emission- 
dominated P-Cygni profiles (Fig. [3J . 



4 GENERAL PICTURE 

The previous discussion of spectroscopic and photometric 
data on SN 1994W (SCL98) led to the conclusion that the 
narrow lines originate in a dense CS envelope. However, at 
present, the spectroscopic and photometric results cannot 
readily be incorporated into any existing model of a SN II 
interacting with a CS environment. We therefore present 
first our qualitative view of what we observe in the case of 
SN 1994W. We emphasize some basic elements of the physi- 
cal picture that, in our opinion, are crucial for understanding 
the phenomenon. 



4.1 The opaque cool dense shell 

The smooth continuum of SN 1994W and its lack of broad 
absorption lines is typical of most SNe Iln at an early phase. 
It is also reminiscen t of the early spectrum of SN 1998S 
jLeonard et al J 120001) . For SN 1998S, the smooth contin- 
uum and the absence of broad absorption lines was ex- 
plained as an effect of an opaque cool dense shell (CDS) 
which formed a t the interfac e of the SN with its circum- 
stellar medium (Chugai 2001). In this situation, an opaque 
CDS is physically equivalent to an expanding photosphere 
with a sharp boundary and without an external, extended 
SN atmosphere. In comparison, the CDS fo rmed in super - 
nova remnants during the radiative stage tPi kelnerl I1954T) 
and in SNe interacting with a moderately dense CS wind 



JChevalier fe Franssoni ri985) are optically thin in the con- 
tinuum. 

The large optical depth of the CDS requires a rela- 
tively large swept-up mass. This may be a natural conse- 
quence of shoc k break-out in SNe II with extended stel- 
lar envelopes JGrasberg et alJ Il97ll : iFalk fe Arnetd Il977l : 
iBlinnikov fe BarfanovlR99^F and/or with extended CS en- 
velopes of unusually high density. In the model of a pre- 
super nova with initial radius ~ 10 15 cm (Fa lk fe Arnetd 
i 1.977?;. a CDS with a mass of ~ 2 Mq forms and re- 
mains opaque for about 70 days. We propose that a similar 
opaque CDS formed in SN 1994W. This suggestion finds 
support in the exceptionally bright peak of the light curve, 
M v » -19.5 mag (SCL98, their Fig. 5). According to the 
theory of SN II light curves, a broad luminous maximum 
with M m — (19—20) mag requires the explosion of a very ex- 
ten ded pre-supernova with an envelope radius of Ro w 10 
cm dGrasberg et alJll97ltlFalk fe Arnetdll977ft . 

4.2 The broad-line region 

The deceleration of the CDS in the extended pre-supernova 
envelope and in the dense CS environme nt is accompa- 
nied by the Rayleigh- Taylor (RT) instability dFalk fe Arnetd 
ll977tlChevaIierll982ft . The growth and subsequent fragmen- 
tation of RT spikes in the dense CDS material results in the 
formation of a narrow layer (AR/R m 0.1 — 0.15) composed 
of dense filaments, sheets and knots embedded in the rar- 
efied hot gas behind the forwar d shock llChevalier fe Blondinl 
Il995t iBlondin fe Ellisonll200ll) . This mixed layer on top of 
the CDS is likely responsible for the broad emission lines ob- 
served in SN 1994W. If this is the case, then the expansion 
velocity of the CDS (v s ) should be the same as the velocity 
of the broad line-emitting gas, i.e., v s ~ 4000 km s _1 . 

A drawback of this picture of the broad-line formation is 
that the velocity distribution of the line-emitting gas peaks 
at the CDS expansion velocity v s , which means the expected 
line profile should be boxy. With the exception of the Na I 
profile on day 57, the broad lines are instead triangular. 

A triangular profile can be produced, however, if the 
spherically symmetric velocity distribution of the line- 
emitting matter is close to dM/dv oc v. Actually, given 
this mass- velocity spectrum and constant emissivity per unit 
mass, and assuming that the velocity field is spherically sym- 
metric, the luminosity distribution over the radial velocity 
u (in units of v s ) is 

dL f 1 ldM 

which is indeed triangular. 

The broad velocity spectrum of the dense gas in the 
mixed layer may naturally arise in two plausible scenarios. 
In the first scenario the CS gas is clumpy. In this case, an 
ensemble of dense CS clouds engulfed by the forward post- 
shock gas could provide a broad velocit y spectrum o f radia- 
tive cloud shocks and cloud fragments lchugailll997l) in the 
range of 1000-4000 km s~ , where the lower limit is the CS 
velocity. This model may be characterized as low-velocity 
clouds in a high- velocity flow. 

Alternatively, dense high-velocity clumps interacting 
with the low-velocity CS gas could produce a similar spec- 
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trum of fragments. This could occur if the forward shock is 
radiative, forming a thin post-shock layer with AR/R < 0.1. 
In this situation, dense RT spikes, formed at the previous 
stage of RT instability of the CDS, can penetrate into the 
pre-shock zone. The interaction of these protrusions with 
the CS gas in the pre-shock zone leads to further fragmen- 
tation and deceleration of fragments down to the velocity of 
the CS gas. As a result, a broad velocity spectrum of dense 
cool fragments emerges in the velocity range of 1000-4000 
km s" 1 . 

For this second possibility to hold, a high CS density 
n is needed to maintain the narrow width of the post-shock 
cooling region, 



32nA' 



(2) 



where m p is the proton mass, and A is the cooling function in 
the post-shock zone. To reach the required ratio l c /R ~ 0.1 
at radius R — 10 15 cm, one needs a CS density of n w 
1.7 x 10 9 cm" 3 assuming A = 2 x 10~ 23 erg s _1 cm 6 and 
v 3 = 4000 km s" 1 . 

Both scenarios for the velocity spectrum leave open the 
question of why the velocity distribution of the line-emitting 
gas in the mixed layer is close to dM/dv oc v. The following 
naive model seems to provide a hint. 

Let us consider the interaction of a high-velocity frag- 
ment with the CS gas. In the rest frame of the initial frag- 
ment, the fragmentation process in the rarefied CS flow 
with the velocity as v s may be thought of as a stripping 
flow with accelerating velocity v\ and in which the "ra- 
dius" a of frag ments progressively decreas es as the velocity 
Vi increases IlKlein. McKee fc Colellalll994f) . Assuming mass 
conservation of fragments in velocity space, Q — dM/dt = 
(dM/dv)(dv/dt) — constant, the mass- velocity spectrum of 
fragments may be expressed as 



dM =Q rdu L y- 

dv\ V dt J 



(3) 



which is essentially determined by the acceleration dvi/dt. 
The latter can be approximately described as acceleration 
due to the drag force 



dvi 2 ,-, x2 

m = na p(l — Vi) . 

dt 



(4) 



where the velocity is in units of v B , m = (4iv/3)a 3 pi is the 
mass of fragment of the density pi, and p is the density in 
the rarefied ambient flow. To determine the acceleration we 
must specify the relation between fragment size and velocity 
v — 1 — vi relative to the CS gas. We assume that the 
fragment size for each new generation of cloudlets is tuned to 
the Kolmogorov turbulent cascade with spectrum v oc a 1 * 13 . 
With this relation, equation J1J leads to dvi/dt oc 1/v. The 
latter combined with equation © leads to a mass-velocity 
spectrum dM/dv oc v. 

A similar analysis carried out for the case of a CS cloud 
in the SN flow produces a quite different mass- velocity spec- 
trum, dM/dv oc (1— v). This spectrum also results in a broad 
line but with zero slope at |u| = 1 and a narrow logarithmic 
peak: dL/du oc (|it| — 1 — In \u\). While this profile resembles 
the observed Ha profile, it is at odds with the triangular 
profile of He I lines. Below (Section 15.41 we argue against a 
large contribution of CS clouds in the formation of Ha. 
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Figure 6. Circumstellar lines on day 79. Overplotted on the ob- 
servational data (dotted line) are the models for homologous ex- 
pansion (thick line) and the constant- velocity case (thin line). 



4.3 The circumstellar envelope 

4-3.1 Expansion kinematics 

The radial velocities of narrow absorption features in SN 
1994W are strikingly persistent between days 18 and 89. 
This implies that the expansion velocity of the CS envelope 
(~ 1000 km s _1 ) has a pre-explosion origin. An alternative 
possibility — acceleration of a slow wind by SN radiation 
(see also SCL98) — seems unlikely, since it would require un- 
realistic fine tuning between the time-dependent luminosity 
and absorption (scattering) coefficient to arrange the persis- 
tence of the velocity of accelerated gas. 

Two extreme options are conceivable for the CS enve- 
lope's pre-explosion kinematics: constant-velocity flow and 
homologous expansion (u oc r). The former might be the 
result of either continuous wind outflow or a succession of 
ejection events with similar ejecta velocities. Homologous 
expansion, on the other hand, might be produced by a sin- 
gle explosive ejection. Th e latter mechanism, though exo tic, 
was in fact proposed by iGrasberg fc Nadvozhinl Jl986l) to 
account for the narrow CS lines in SN 1983K. 



4-3.2 Envelope size 

The light curve of SN 1994W is characterized by a broad 
maximum around day 30 and a subsequent plateau, ter- 
minated by a sudden drop in luminosity at id ~ 110 d 
(SCL98, their Fig. 5). Assuming a CDS expansion velocity of 
v s — 4000 km s _1 , we find that during this period the CDS 
sweeps up the CS gas within the radius Rd = v B td ~ 4 x 10 15 
cm. Although this estimate refers to t = 110 d, we adopt it 
as a rough estimate for earlier epochs as well. 

We now argue that this radius should coincide with the 
outer radius of the dense CS envelope, R cs ~ Rd- If in- 
stead R CB > Rd, then the light curve of SN 1994W should 
not have shown such a steep drop, since the CS interaction 
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would have augmented the luminosity at the luminosity de- 
cay phase. Equally, if _R CS < R<\, the narrow lines would have 
disappeared at epoch t < t c \, which was not observed. These 
arguments, therefore, provide strong evidence that the ex- 
tent of the dense CS envelope is R CB ~ Rd ~ 4 x 10 15 cm. Be- 
yond this radius, the CS density presumably drops steeply, 
as indicated both by the broad-band light drop after day 
100 and the fast decline of the Ha luminosity (Fig. . At 
larger radii, X-ray observations suggest t hat other den sity 
enhancements may be present (Section 171 ISchlegellll999|) . 



4-3.3 Electron-scattering wings and circumstellar density 

We attribute the extended smooth wings observed in Ha and 
H/3 between days 18 and 89 to the effect of Thomson scatter- 
ing on thermal elect rons participa ting in the bulk expansion 
of the CS envelope JChugaill200l . Based on previous com- 
putations of Thomson scattering effects in SN 1998S, and 
the strength of the Ha wings in SN 1994W on day 57 (when 
blending with He I A6678 was negligible), we believe that 
the optical depth of the CS envelope to Thomson scattering 
(tt) at this stage must have been close to unity. 

Adopting a photospheric radius R p » v B t = 2 x 10 cm 
on day 57 (where v s « 4000 km s _1 ) and an outer radius of 
the CS envelope R CB — 4 x 10 15 cm, we obtain an estimate 
of the average electron concentration in the CS envelope 
assuming tt ~ 1: 



0-t(Rcs — Rp) 



10 9 cm" 3 . 



(5) 



This value, taken together with the forward shock- wave ve- 
locity (v B h ~ 3000 km s _1 ), after correction for the CS ve- 
locity, provides an estimate of the density of the CDS mate- 
rial in the broad-line region (n c d s ). The pressure equilibrium 
condition, with sound speed in the CDS matter of c s « 15 
km s _1 , results in n c d B ~ n(u s h/c s ) 2 ~ 4 x 10 13 cm -3 . This 
tremendous density provides a natural reason for the strong 
collisional thermalization of the broad component indicated 
by the observed inverse Balmer decrement (Section 13.3. II . 



4.3.4 Circumstellar density from narrow lines 

The strength of narrow subordinate lines indicates a high 
density in the CS envelope, as noted by SCL98. The sim- 
plest density estimate may be taken from the condition that 
the optical depth in a line of Fe II, for example A5018, is on 
the order of unity. This follows from the relative intensity 
of the absorption component of this line (~ 0.5). Assuming 
that level populations of Fe 11 obey the Boltzmann distribu- 
tion for T = 10 4 K, one obtains for solar Fe/H a hydrogen 
concentration of n ~ 3 x Vs v%r^ b cm -3 . Here ris is the 
linear scale of the line-forming zone in units of 10 15 cm, vs 
is the velocity dispersion in units of 10 s cm s , and /a is 
the ionization fraction of Fe II. Given the scale of the CS en- 
velope (~4x 10 15 cm), we thus have a lower limit n > 10 7 
cm -3 . This lower limit should be considered as a revised 
version of the estimate reported in SCL98. 

A somewhat more elaborate estimate of the density 
may be obtained using several different lines and assuming 
a quasi-local thermodynamic equilibrium (quasi-LTE) ap- 
proximation, i.e., the Saha-Boltzmann equations corrected 



for both geometrical dilution W and dilution of the pho- 
tospheric brightness f. The latter defines the photospheric 
brightness through the black-body brightness I v = £B„(T). 
We adopt a simple density distribution in the CS envelope: 
a plateau with a steep outer drop 



n = n /[l + {r/R^f 



(6) 



where the cutoff radius i?k ~ Res = 4 x 10 cm. We 
consider two cases for the CS kinematics: free expansion 
(v = hr), with outer velocity at r — 6 x 10 l0 cm equal to 
Vb = 1100 km s , and constant-velocity flow with v = 1100 
km s . The photospheric radius on day 79 is estimated as 



V s t 



2.7 x 10 b cm, where v s = 4000 km s , while 



the temperature is taken to be T p — 7200 K according to 
the value found from the black-body fit to the continuum 
(Section 13.2.31 . The parameter £ in this case is found to be 
0.26. 

In Fig.|B]we show calculated profiles of Fe 11 A5018, Si 11 
A6148, the Na I Di j2 doublet and the O I 7773 A triplet, for 
solar abundance and hydrogen concentration no = 1.5 x 10 9 
cm -3 . In the case of Na I we have for the sake of simplicity 
ignored the doublet structure, which explains why our model 
line is narrower than the observed one. This density seems 
to be the optimal one: a factor of 1.5 lower density results in 
too weak O I and Na I lines, while a factor 1.5 higher density 
makes all the lines too strong. 

Although neither kinematic model produces an excel- 
lent fit, it is clear that free expansion kinematics better pre- 
dicts the positions of the absorption minima and, therefore, 
is preferred compared to the constant-velocity case. The dif- 
ferences between the modelled and observed line profiles in 
the free-expansion case may be related to the omission of 
Thomson scattering and the simplicity of our quasi-LTE 
model. 



4.4 Overview of the qualitative model 

We summarize the main results of the qualitative analysis 
in a cartoon (Fig.0. It shows the basic structural elements 
of SN 1994W which we believe are responsible for the for- 
mation of the optical spectrum when the supernova's lumi- 
nosity is high. The ejecta expand with a velocity of ~ 4000 
km s _1 into an extended (~4x 10 14 cm) CS envelope. The 
characteristic density of the CS envelope is ~ 10 9 cm -3 . The 
SN ejecta are enshrouded by the opaque CDS within which 
the photosphere resides during most of the luminous phase. 
The opaque CDS precludes the formation of absorption lines 
from the ejecta. The CDS presumably forms primarily dur- 
ing shock break-out and is subsequently maintained by both 
CS material swept up by the radiative forward shock and SN 
material swept up by the reverse shock. These shocks are not 
shown in Fig. Q but they are presumably located at the dis- 
tances AR < 0.1R away from the outer and inner edges of 
the CDS. 

Attached to the CDS is a layer populated by dense 
(n ~ 4 x 10 13 cm -3 ) fragments supplied by the RT insta- 
bility of the dense CDS and, possibly, by radiative shocks 
in CS clouds. This inhomogeneous layer of dense material is 
the primary site of the broad emission lines with the charac- 
teristic velocity of the line-emitting gas, 4000 km s _1 (e.g., 
in He 1). 
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Figure 7. A sketch of SN 1994W at the epoch around day 30. We 
show the principal structural elements involved in the formation 
of the spectrum. The SN ejecta are bounded by an opaque cool 
dense shell (CDS), which is responsible for the continuum radi- 
ation. The broad-line region is a narrow mixing layer attached 
to the CDS, composed of Rayleigh- Taylor fragments of the CDS 
matter and possibly of shocked CS clouds. The SN ejecta expand 
into a dense CS envelope with Thomson optical depth on the or- 
der of unity. The CS envelope is responsible for both narrow lines 
and the extended Thomson wings seen in Ho. 

The CS envelope, which in turn expands with a veloc- 
ity of ~ 1000 km s" 1 , is responsible for the narrow lines. 
Thomson scattering in the CS envelope results in the emer- 
gence of broad emission-line wings, which are most apparent 
in the strong lines such as Ha and H/3. Because of the large 
expansion velocity, the red wing is stronger than the blue 
wing. 



5 A MODEL FOR Ha 

In order to confirm and refine the above qualitative picture 
of the formation of the SN 1994W spectrum, we have mod- 
elled the Ha line profile. In many respects our a pproach 
repeats the analysis of SN 1998S bv IChueail (1200 ll) . which 
involved calculations of the CDS dynamics and Monte Carlo 
computations of the emergent lines formed outside the CDS. 
Here we use an updated version of this model, which now in- 
cludes calculations of CS hydrogen ionization, electron tem- 
perature and Ha emissivity. In the previous model both of 
the latter were set rather arbitrarily. 

5.1 The model 

The CDS dynamics are calculated nu merically in the thin- 
shell approximation (IChevalierl 1 19821) . We assume that the 
SN initially expands homologously (v ex r) with density dis- 
tribution p — po/[l + (v/Vk) 9 ], where the parameters of the 
distribution are defined by the ejecta mass M and kinetic en- 
ergy E. At radius Rq = 1000 Rq (a rather arbitrary value) 



the SN ejecta begin to interact with the extended stellar 
envelope. The density of the extended stellar envelope is ad- 
justed in our model by the requirement that a strong initial 
deceleration of the swept-up thin shell must result in an 
expansion velocity of ~ 4000 km s _1 on day 30 when the 
photospheric radius reaches ~ 10 15 cm. The CS envelope is 
attached to the stellar envelope at r m 10 cm. A constant 
pre-shock velocity (1000 km s _1 ) is assumed for the CS gas. 
This is a reasonable approximation for the homologous ex- 
pansion, bearing in mind additional radiative acceleration. 
For the SN mass and energy we adopt M = 8 Mq and 
E — 10 51 erg — both values are cl ose to the parameter s 
of the light-curve models B and C of lFalk fc Arnetd <1977t) . 
We find that the dynamical model is nearly the same for SN 
mass in the range 7—9 Mq. The evolution of the CDS radius 
and velocity is shown in Fig.|H|for the density profile which is 
also consistent with the Ha model. In fact, the optimal den- 
sity is a compromise between the requirements imposed by 
the Ha profile evolution, the CDS velocity, and the bolomet- 
ric luminosity at the final stage of the light-curve plateau. 
The density of the CS envelope is nearly flat (p oc r~ 0,4 ) in 
the range 10 15 < r < 4 x 10 15 cm. 

The ionization of the CS envelope is produced by both 
X-ray emission from the radiative forward shock wave and 
by photospheric radiation. The latter primarily operates via 
hydrogen photoionization from the second level. The X-ray 
luminosity of the forward shock wave is calculated as 

L x = 2nr 2 p(v s ~ it ps ) 3 , (7) 

where p is the pre-shock density, while the pre-shock ve- 
locity u ps is the superposition of the pre-explosion velocity 
and the accelerated velocity term. Half of this luminosity is 
directed toward the CDS where it is reprocessed into con- 
tinuum radiation. This component is not treated explicitly 
in the model. The other half of the X-ray luminosity of 
the forward shock is emitted outward and is partially ab- 
sorbed by the CS gas. The X-ray spectrum of the forward 
shock is approximated as Fe = CE~ q exp (— E/T s ), where 
T s is the temperature of the forward shock defined by the 
shock velocity v s — u ps . The spectral power index q — 0.5 
roughly describes both the Gaunt factor an d the contribu- 
tion of X-ray lines in the low-energy band jTerlevich et alJ 
1992). The soft X-rays from the reverse shock at the epoch 
und er consideration are fully absorbed by t he CDS (cf. Fig. 
9 of lFransson. Lundavist fc Chevalierlll996l) and by the SN 
ejecta, and are reprocessed into the optical. 

The energy of the X-rays absorbed in the CS envelope 
is shared between heating, ionization and excitation. Ex- 
citation and ionization by the photospheric radiation are 
also taken into account. The electron temperature T c in the 
CS envelope is determined from the energy balance between 
heating from ionization by X-rays and photospheric radi- 
ation and cooling due to hydrogen. This approximation is 
quite reasonable in the temperature range of (1 — 2) x 10 4 
K. The hydrogen atom is treated in the two- level-plus- 
continuum approximation. The degree of ionization we ob- 
tain for hydrogen in the CS envelope lies typically in the 
range 0.8 — 0.95. The Ha emissivity is due to recombination 
(Case B) and collisional excitation. Collisional de-excitation 
of Ha is included in the computation of the net emissivity. 
To allow for uncertainties in the flux, model, distance and 
reddening we use a fitting parameter C n for narrow Ha. This 
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Figure 8. The CS density {top panel) and evolution of the radius 
and velocity of the cool dense shell {lower panels). 



parameter is a multiplicative factor in the Ha emissivity in 
the CS envelope. 

The broad-line region (Fig. is described by an (on 
average) homogeneous layer, R s < r < R\, where Ri/Rs = 
1.15 and the distribution of the line luminosity dL/dv oc v 
for vi < v < v s , where vi is the minimal velocity adopted to 
be equal to the pre-shock CS gas velocity. We express the 
average emissivity of the broad component in terms of the 
pre-shock emissivity, multiplied by a fitting factor Cb, which 
is restricted by the requirement that the Ha intensity from 
the broad-line region cannot exceed the black-body value for 
the local electron temperature of the CS gas in the pre-shock 
zone. 

We assume free expansion pre-explosion kinematics 
(u = hr), though we have also explored the constant- velocity 
case. The acceleration by the SN radiation is described as 
an additional CS velocity term u a = uo{R s /r) 2 , where uq is 
a free parameter. 



5.2 Thomson scattering and the broad component 

In the Ha line profile, the broad emission component and 
electron-scattering wings overlap and at first sight cannot 
be disentangled unambiguously. Nevertheless, we found that 
there is not much freedom in the decomposition procedure. 
In Fig.|5]we show a model narrow-line component on day 31, 
with the broad component turned off. To emphasize the dif- 
ferent behaviours of the broad component and the electron- 
scattering wings, we show in Fig. I1UI models with different 
Thomson optical depth (0, 1 and 2.8) but with otherwise 
similar parameters. For optical depth zero, the modelled 
broad component is strongly skewed toward the blue be- 
cause of occultation by the photosphere. This asymmetry 
decreases as tt increases and Thomson scattering of the 
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Figure 9. A typical model of He? on day 31 {thick line), showing 
the contribution of the narrow component ( thin line) . 
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Figure 10. The same model as in Fig.EObut for different Thom- 
son optical depths. Lines with increasing thickness show cases 
it = 0, 1 and 2.8, respectively. 

line emission in the expanding CS envelope produces an in- 
creasingly strong red wing. These different behaviours of the 
broad line and electron-scattering wings permit us to disen- 
tangle the contribution of the two line components in Ha. 
In passing, we note that the continuum level also depends 
on tt, since the Thomson optical depth affects the escape 
probability of photons emitted by the photosphere. 

5.3 Ha and the expansion law 

The sensitivity of the model line profile to the kinematics 
of the CS envelope is demonstrated in Fig. 1111 The plot 
shows Ha profiles on day 31 computed for similar model 
parameters — all that differs is the velocity distribution. 
For tt = 2.8, homologous expansion without post-explosion 
radiative acceleration fits the observations quite well. The 
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Figure 11. Ha profile on day 31 for different velocity distri- 
butions (shown in the inset in each panel). Overplotted on each 
observed profile (thin line) are models for CS envelope kinematics 
displayed in insets. 

constant-velocity case is much less successful: it produces 
an emission component which is too broad, and absorption 
which is too shallow (Tig. 111)1. This result, taken together 
with the results of the narrow CS line modelling (Fig.|SJ, ar- 
gues against the constant-velocity case. We checked a model 
with constant velocity and pre-shock acceleration and found, 
unsurprisingly, that this gives even worse agreement. The 
homologous model, on the other hand, works just as well 
when combined with post-explosion acceleration. In the bot- 
tom panel of Fig. 1111 we show such a model profile char- 
acterized by an amplitude of accelerated velocity term of 
u a = 400 km s" 1 . 

Although homologous expansion kinematics is certainly 
preferred, this model shows a rather deep absorption compo- 
nent. This contradiction is not related to the finite resolution 
since we have smoothed the model profile by a Gaussian with 
appropriate width. An explanation for the large absorption 
strength in the model might be hidden in a possible devi- 
ation of the kinematics from homologous expansion, some 
dumpiness of the CS envelope, or lower hydrogen excita- 
tion in the outer layers of the CS envelope than the model 
predicts. 

5.4 Modelling the Ha evolution 

We have calculated Ha for three epochs (31, 57 and 89 days 
past explosion) using the density distribution of the CS en- 
velope (r > 10 15 cm), and the CDS radius and velocity 
as shown in Fig. |S| The density parameter w = Aur 2 p at 
r = 10 15 cm on day 31 is w = 7.5 x 10 16 g cm' 1 . Free ex- 
pansion kinematics is assumed with a boundary velocity of 
1100 km s _1 at r = 5.4 x 10 15 cm on day 31. This implies 
an age for the CS envelope of t cs = 1.5 yr when the SN 
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Figure 12. Ha at different epochs. Model profiles (thick line) are 
overlaid on the observations (thin line). 



explodes. For other values of t cs counted from the moment 
of envelope ejection, we rescale according to r a t cs and 
P « 1/tls. 

The computed profiles are shown in Fig.H"2l In Tabic H 
we list the corresponding parameters: photospheric (CDS) 
radius, photospheric temperature, brightness dilution, CDS 
velocity, post-explosion velocity increase, and emissivity fit- 
ting factors for narrow and broad components. The last col- 
umn displays the calculated Thomson optical depth outside 
the CDS. On days 31 and 57 we use the approximation of an 
opaque photosphere. On day 89 we found, however, that this 
approximation predicts too weak a red part of the profile. 
The model agrees better with the observations if we suggest 
that at this late phase the CDS is semi-transparent. To de- 
scribe this effect in a simple way we allow a photon striking 
the photosphere to cross it and escape into the CS medium 
with a finite probability, which we found should be close to 
p = 0.3. This assumption is qualitatively consistent with the 
fact that the photospheric brightness on day 89 is diluted 
(£ = 0.2) for the optimal continuum temperature. The mod- 
els for all three epochs demonstrate satisfactory fits to the 
data. On day 31, the model is consistent with the presence 
of broad He I A6678 in the red wing of Ha, while on day 57 
the model fit in the red wing does not leave any room for the 
He I line. This behaviour is consistent with the weakening 
of other He I lines at this time. 

We calculated the Ha component from the CS envelope 
in a unified way with a minimum number of free parameters. 
For this reason, the fact that the tuning parameter of the 
narrow- line intensity (C n ) is close to unity (Table means 
that within the uncertainties (~ 20 per cent) the electron 
distribution in the CS envelope recovered from the Thomson 
scattering effects is also consistent with the luminosity of the 
narrow Ha component. Because the line luminosity depends 
on the filling factor as L oc 1/f while the Thomson optical 
depth does not, we conclude that the filling factor in the CS 
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envelope is close to unity in at least most of the line-forming 
zone. 

This last observation is of crucial importance for distin- 
guishing between the two scenarios for the formation of the 
broad- line region (Section l4,21 , and indicates that CS clouds 
do not play a significant role in this process. 

We have studied cases with different density power-law 
indices s (defined via p cc r s ) in the flat part of the CS 
envelope. In the case s = 0, the tuning parameter (C n ) must 
systematically decrease by a factor of 1.25 to fit the profile. 
For s — — 1, C n should systematically increase with time by 
a similar factor. Bearing in mind the uncertainties in the 
modelling and fluxing, we do not rule out that the density 
gradient lies in the range — 1 ^ s ^ 0. 



5.5 Implications for the origin of the 
circumstellar envelope 

The model for the Ha evolution presented above for the 
CS envelope outside r ~ 10 15 cm argues for a total mass 
of M cs M 0.4 Mq and kinetic energy of E CB « 2 x 10 48 
erg. Combined with the derived kinematic age of the CS 
envelope t cs m 1.5 yr, the average kinetic luminosity of the 
mass-loss mechanism responsible for the ejection of the CS 
envelope then becomes L m E cs /t cs r;4x 10 40 erg s"\ and 
the average mass-loss rate is ~ 0.3 Mq yr -1 . 

The kinetic luminosity we derive here is enormous com- 
pared to stellar values. It exceeds by two orders of mag- 
nitude the radiative luminosity of a pre-supernova with a 
main-sequence mass of ^ 20 Mq. This certainly rules out 
a superwind as the mass-loss mechanism responsible for the 
CS envelope around SN 1994W. The CS envelope around 
SN 1994W must have been born as a result of a rather vio- 
lent mass ejection initiated by some energetic explosive event 
in the stellar interior (E ^ 2 x 10 48 erg) at ~ 1.5 yr prior to 
the SN outburst. 



6 LIGHT-CURVE MODELLING 

To model the broad-band photometric light curves, we 
use the multi-energy group radiatio n hydrodynamics code 
STELLA ijBlinnikov et al.lll99Sl . l200oT) . In the current work, 
STELLA solves time-dependent equations for the angular mo- 
ments of intensity averaged over fixed frequency bands, us- 
ing 200 zones for the Lagrangian coordinate and up to 100 
frequency bins with variable Eddington factors. The transfer 
of gamma rays from radioactive decay is calculated using a 
one-group approximation for the non-local deposition of the 
energy of radioactive nuclei. Except for the latest phases, 
the gamma-ray deposition is not important for SN 1994W. 
In the equation of state, LTE ionizations and recombina- 
tions are taken into account, but radiation is not assumed 
to be in equilibrium with matter. The effect of line opac- 
ity is treate d as an expansion o pacity according to the pre- 
scription of lEastman fc Pintol lll993t) . In comparison with 
previous work with STELLA, here we use an extended spec- 
tral line list, t he sa me as in the discussion of SN 1998S by 
IChueai et all (I2002I) . 
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Figure 13. Composition as a function of interior mass, M r , for 
the most abundant elements in the pre-supernova model sn94w58. 
The mass cut is at M c = 1.41 Mq. 



6.1 Pre-supernova models 

SCL98 used anal ytical fits jLitvinova fe Nadvozhinl Il983l . 
119851: lPopovl ll993') to estimate the pre-SN radius, mass and 
explosion energy in SN 1994W. However, those fits are only 
valid for typical SNe II-P, and from the discussion above we 
know that SN 1994W was instead dominated by circumstel- 
lar interaction. To determine the pre-supernova parameters 
for such an unusual event, we therefore need to do detailed 
numerical modelling. This can be done only for a subset of 
the pre-supernova parameters, so we used a cut-and-trial ap- 
proach, guessing the initial conditions and computing light 
curves in order to meet the constraints posed by the photom- 
etry and the spectral analysis of Section^ In the absence 
of an evolutionary model, we built a sequence of several 
dozens of non-evolutionary models, converging finally on a 
reasonably good fit to observations. Our initial models are 
constructed in hydrostatic equilibrium for the bulk mass, 
while the outer layers mimic the structure of the circumstel- 
lar envelope. 

For a given mass M and radius Ro , we obtain a model in 
mechanical equilibrium assu ming a power-law dependence 
of temperature on density jNadvozhin fc Razinkoval Il986l : 
iBlinnikov fc Bartunovlll993h : 

T oc p a . (8) 

The hydrostatic configuration thus obtained would be close 
to a polytrope of index 1/a if it were chemically homoge- 
neous and fully ionized. The difference from a polytropic 
model arises due to recombination of ions in the outermost 
layers and non- uniform composition (Fig. I13t . 

In the centre of this configuration, at the mass cut of the 
collapsing core, we assume a "point-like" gravitating hard 
core (with numerical radius R c = O.IRq — much larger 
than a real core, but much smaller than the radii of mesh 
zones involved in our hydrodynamic simulations). The den- 
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Table 4. Parameters of the Ha evolution models 
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Figure 14. Density as a function of the radius r in the pre- 
supernova models sn94w43 and sn94w58. 



sity structure found in this way is shown in Figs. 1141 and 

The parameters of the models are given in Table 16.11 
The first column is the model label, which is followed by the 
mass of SN ejecta, pre-supernova stellar radius, Ni mass, 
power-law index a, density of the CS envelope (pis) at the 
radius r = 10 15 cm, power-law index (p) of the CS envelope 
density distribution p oc r~ p , the outer radius of the CS 
envelope, and the kinetic energy at infinity in units of foe 
(1 foe = 10 51 erg). 

Each model was exploded by the deposition of heat en- 
ergy in a layer of mass ~ 0.06 Mq outside of 1.41 Mq. Since 
STELLA does not include nuclear burning, preservation of the 
same mixed composition in the ejecta is assured. 

We explored parameter space for the mass (M) and en- 
ergy (E) of the SN ejecta and found acceptable fits to the 
data for masses in the range 6-15 Mq and with the en- 
ergy given by the ratio E/M « 0.15-0.2 foe Afg 1 . In the 
following sections and in Table 16.11 we present three rep- 
resentative models which allow us to demonstrate different 
aspects of the formation of the light curve. Of these three, 
model sn94w58 provides our best fit to SN 1994W, model 
sn94w43 shows that a normal SN II-P light curve is not 
appropriate for SN 1994W, and model sn94w68 illustrates 
what the supernova light curve might have looked like with- 
out CS interaction. All three models have E = 1.5 x 10 51 
ergs. This energy is typical for good light-curve models. The 




5 10 15 

M r /M 



Figure 15. Density as a function of the interior mass (M r ) in the 
pre-supernova models sn94w43 and sn94w58. The mass cut is at 
M c = 1.41 Mq. 

asymptotic kinetic energy of the ejecta is somewhat lower 
and is given in the Table ISTl 

6.2 Hydrodynamics and shock propagation 

The modelled light curves are dominated by the diffusion 
of the trapped radiation generated during the shock wave's 
propagation in the extended stellar atmosphere and sub- 
sequent emission of a radiative shock propagating in the 
dense circumstellar medium. A dense shell is formed which 
is found in non-equilibrium radiation hydrodynamic mod- 
elling but often missed in equilibrium diffusion modelling. 
We identify this dense shell as the CDS discussed in Sec- 
tion 14.11 The Lagrangian code STELLA does a good job of 
resolving the very fine structure of the opaque shell, which 
contains approximately one solar mass of the material (Fig. 

6.3 Light curves 

In Figs. I17l and ll8l we show the changes in the colour tem- 
perature, T c , of the best black-body fit to the flux. We 
compare this to the effective temperature, T e g, defined by 
the luminosity and the radius of last sca ttering R through 
L = 4-KaT* B R 2 (see lBlinnikov et aLlllifflsl for details of find- 
ing R and from that T c g). Our multi-group radiative trans- 
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Table 5. Parameters of hydrodynamical models 



Model 


M ci a 


Ro 


Mm 


a 


Pl5 


v Rv> 


Skin'' 




(Mq) 


(10 4 R e ) 


(Me) 




(10~ 15 g cm" 3 ) 


(io 4 Re) 


(10 51 erg) 


sn94w43 


12 


8 


0.015 


0.319 


5 


2 12 


1.07 


sn94w58 


7 


2 


0.015 


0.31 


12 


1 6.6 


1.25 


sn94w64 


7 


2 


0.015 


0.31 







1.33 



a pre-supernova mass = M e ; + 1.41 Mq 
^ kinetic energy at infinity 
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Figure 16. Density as a function of radius (r) in two models with 
different amounts of artificial smearing at day 30. The height of 
the density spike with small smearing (solid line) reaches the 
theoretical estimate for an isothermal shock wave. 



fer with hydrodynamics obtains this temperature in a self- 
consistent way, and no additional estimates of the thermal- 
iz ation depth are needed (in contrast to the one-group model 
of lEnsman fc Burrowslll99ol for example). The large differ- 
ence between colour and effective temperatures is partly due 
to a geometric effect (the radius of last scattering is greater 
than the effective radius of photon creation) and partly due 
to the blanketing effect of scattering (the average energy of 
the photons is high er than that corresponding to the value of 

ISchusterll9rja) . The modelled T c changes in a way sim- 
ilar to the values derived from black-body fits to the spectra 
(Section 13.2.31 . 

The light curves are shown in Figs. E3 and|2T] The 
model sn94w43 produces a bright light curve, but the sec- 
ond half of its plateau is dominated by diffusion, not by the 
shock, so the behaviour of colours is similar to that of a 
typical SN II-P. 

The density of the more powerful wind in the model 
sn94w58 follows the law p oc r . This model produces a 
much better fit to the colours at the late plateau phase. The 
flat CS density distribution is consistent with the modelling 
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15 



10 




50 100 150 

t, days since 1994 July 14.0 

Figure 17. Colour and effective temperatures for the model 
sn94w43. Realistic, scattering-dominated opacity has been as- 
sumed. The solid line shows the temperature of the best black- 
body fit to the flux (colour temperature). The dashed line shows 
the effective temperature defined by the luminosity and the radius 
of last scattering. 



of Ha evolution (Section 15.41 . Model sn94w64, which is ex- 
actly the same as model sn94w58 in the bulk mass, but does 
not have the powerful wind, is much less luminous at the late 
plateau phase. This shows the importance of CS interaction 
for the luminosity. 

The time-dependence of the radius R s and velocity V s 
of the cool dense shell in model sn94w58 is shown in Fig. 
1221 The two nearly indistinguishable lines in the V s plot are 
obtained by taking V s = dR s /dt numerically, and by taking 
the mass-averaged speed of matter inside the shell. A com- 
parison with Fig.|S|shows good agreement with the spectral 
model, although V s is somewhat higher in the hydrodynam- 
ical model. The difference is at most ~ 10 per cent. 

To summarize, the hydrodynamical modelling of the 
light curve of SN 1994W suggests that the optimal light 
curve is produced in a model which contains a dense CS en- 
velope at r > 10 15 cm with a relatively flat density distribu- 
tion and an outer radius of ~ 4.5 x 10 15 cm. This conclusion 
is fully consistent with the Ha modelling. 
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Figure 19. Light curves for the model sn94w43. The photometric 
points are taken from SCL98 and references therein. 



7 DISCUSSION AND CONCLUSIONS 

We have presented and analysed spectra and light curves of 
SN 1994W, one of the best-observed SNe Iln. During the 
first three months the spectrum clearly shows the presence 
of three line components: (i) narrow P-Cygni lines, (ii) in- 
trinsically broad lines and (iii) extended smooth wings in Ha 
and H/3. We attribute these components, respectively, to (i) 
a dense CS envelope, (ii) shocked, cool, dense gas confined 
in a narrow layer on top of the photosphere, and (iii) the 
effects of Thomson scattering in the CS envelope. 

Our line profile analysis and hydrodynamical light- 
curve modelling have led us to a coherent picture. SN 1994W 
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Figure 21. Same as in Fig. 1191 but for model sn94w64. 

appears to have been the result of explosion of an extended 
pre-supernova (~ 10 15 cm) embedded in an extended CS 
envelope [~ (4 — 5) x 10 15 cm] with Thomson optical depth 
~ 2.8 one month after the explosion. 

Although we do not rule out the presence of a wind from 
a normal red supergiant pre-supernova outside the dense CS 
envelope, such a wind component cannot have been very 
dense, considering the low luminosity of Ha on day 203 
(L ~ 10 38 erg s _1 ). This luminosity is two orders of mag- 
nit ude lower than for SN 1979 C (~ 10 40 erg s _1 ) at ~ 1 
yr llChevalier fc Fra nsson 1985). Attributing this difference 
to the difference in the interaction luminosity (L tx Aitpvj ), 
we can estimate the density parameter of the outer wind in 
SN 1994W. The expansion velocity of the CDS in SN 1994W, 
~ 4000 km s _1 , is half that of SN 1979C (~ 8000 kms" 1 ). 
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Days 

Figure 22. The evolution of the radius and velocity of the cool 
dense shell in hydrodynamical model sn94w58. Compare with 
Fig.! 



Given the SN 1979C wind parameter, w « 10 16 g cm - , 
the outer wind in SN 1994W should have w ~ 10 15 g cm -1 
to account for the difference in the Ha luminosity. This is 
at least ten times lower than in SN 1979C, and a factor of 

Lundavist & Fransson 
1988). In this respect, the detection bv ISchlerelTi^99h of 
X-rays from SN 1994W at t = 1180 days with luminosity 
L « 8 x 10 39 erg s _1 appears odd when compared with the 
expected X-ray luminosity ( < 10 38 erg s -1 ) at that epoch 
for a case of SN/wind interaction with w » 10 15 g cm -1 . 
To resolve this problem, one would need the wind density 
to be at least of factor of ten higher at ~ 4 x 10 16 cm, i.e., 
exceeding that of SN 1979C. This could argue for a possible 
multi-shell ejection scenario prior to the explosion. 

The recovered kinematics, density and linear scale of the 
CS matter around SN 1994W imply a kinematic age for the 
CS envelope of t cs ps 1.5 yr, mass M CB ~ 0.4 Mq and kinetic 
energy _B CS ~ 2 x 10 48 erg. The enormous average mass-loss 
rate (M cs /t cs ~ 0.3 Mq yr -1 ) and equally enormous kinetic 
luminosity (E CB /t CB ~ 4 x 10 40 erg s _1 ) of the mass-ejection 
mechanism strongly suggest that the CS envelope has been 
lost as a result of an explosive event which occurred ~ 1.5 
yr prior to the SN outburst. 

In their model fo r the narrow lines in SN 1983K, 
iGrasberg fc Nadvozhinl il986l) invoked the violent ejection 
of the hydrogen envelope 1-2 months before the SN explo- 
sion with a velocity on the order of ~ 2500 km s _1 and an 
energy (1 — 2) x 10 49 erg. Dramatically high mass loss shortly 
be fore explosion has also been derived for the SN Iln 1995G 
by IChueai fc Danzigerl )().'!) , on the basis of an analysis 
similar to the one we have presented here. 

The mechanism behind such violent ejections might 
be associated with a nuclear flash in the degenerate core 
of the pre-supernova. Th is conjecture, also presented by 
IChueai fc Danzigerl (120031). is prompted by the prediction 
of IWeaver fc Wooslevl (Il979l) that the O/Ne/Mg core of an 



~ 11 Mq star may produce strong Ne flashes several years 
prior to the SN explosion and that the strongest flash could 
eject most of the hydrogen envelope with velocities of ~ 100 
km s _1 . Although such fl ashes do not occur in more recen t 
models with finer zoning dWooslev. Heeer fc WeaTeill2002h . 
we believe that the effect warrants further investigation, es- 
pecially given the highly complicated nuclear burning regime 
and hydrodynamics of degenerate O/Ne/Mg cores. 

Remarkably, the low 56 Ni mass (< 0.015 Mq) estimated 
from the tail i?-band luminosity of SN 1994W (SCL98) 
seems to be in accord with the suggested mass of such a 
progenitor, if at the lower end of our estimated mass range. 
Such stars are expected to eject only small amounts of 56 Ni 
jMavle fc Wilsonlll988l) . Our models are, however, not very 
sensitive to the zero-age main-sequence mass of the progen- 
itor, and, as SCL98 point out, the low nickel content could 
indicate either a low-mass progenitor which may have lost a 
few solar masses in shell and wind ejections, or a more mas- 
sive star that lost several solar masses prior to the ejection of 
the CS shell. Deep X-ray searches for more extended CS ma- 
terial could help solve this problem, as should more detailed 
models for stellar evolution just prior to core collapse. 

How common are events like SN 1994W? Of the SNe Iln 
that have been observed so far, it is possible that a close 
counterpart has been seen but not recognised as such due 
to sparser temporal coverage or greater distance. Never- 
theless, several objects have shown narrow hydrogen and 
Fe II absorption lines with velocitie s of about 10 00 km 
Thi s family includes SN e 1987B dSchleeel et alJ 



1994aj jBenetti et alJ Il99£ft. 1994ak jFilippenkcl 



1996), 



1997T 

1995G iPastorello et alll2002l). 1996L feenetti et al.lll99' 
and 1999el iDi Carlo et all2002l) . We suggest that in at least 
some of these SNe, a CS envelope was ejected in a violent, 
explosi ve manner, as we beli e ve wa s the case for SN 1994W, 
and as IChugai fc Danzigerl B2003T) have suggested for SN 
1995G. Further data and analyses are needed, however, be- 
fore it can be demonstrated that the energy of the ejected 
CS envelope is comparable for these events. 

In any case, we suggest that at least part of the variety 
of SNe Iln has been accounted for. Some, like SN 1988Z and 
SN 1998S, have a dense, slow CS envelope formed by a slow 
superwind. Others, like SN 1994W, seem to have ejected a 
CS envelope in a violent event a few years before the SN 
explosion. 
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